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We revisit the question of the imposition of initial data representing astrophysical gravitational 
perturbations of black holes. We study their dynamics for the case of nonrotating black holes by 
numerically evolving the Teukolsky equation in the time domain. In order to express the Teukolsky 
function \P explicitly in terms of hypersurface quantities, we relate it to the Moncrief waveform 
<j)M through a Chandrasekhar transformation in the case of a nonrotating black hole. This relation 
between 9 and 4>m holds for any constant time hypersurface and allows us to compare the com- 
putation of the evolution of Schwarzschild perturbations by the Teukolsky and by the Zerilli and 
Regge- Wheeler equations. We explicitly perform this comparison for the Misner initial data in the 
close limit approach. We evolve numerically both, the Teukolsky (with the recent code of Ref. 
and the Zerilli equations, finding complete agreement in resulting waveforms within numerical error. 
The consistency of these results further supports the correctness of the numerical code for evolving 
the Teukolsky equation as well as the analytic expressions for ^ in terms only of the three-metric 
and the extrinsic curvature. 



I. INTRODUCTION AND OVERVIEW 



Binary black holes provide, in principle, one of the strongest sources of the gravitational radiation to be observed 
by the detectors currently under construction. Black holes have proved to be very elusive to standard astronomical 
methods of detection. This is a consequence of the fact that astrophysical black holes, being one of the simplest 
objects in nature, interact only gravitationally with the rest of the universe. Thus, gravitational wave observatories 
are particularly suitable to detect them. The confirmation of black hole existence in nature can be considered to be 
as great a discovery as the confirmation of the existence of gravitational radiation. 

This exciting new experimental-observational situation has taken theory a bit by surprise. In the early seventies 
a great deal of effort has been devoted to study gravitational perturbations by means of the Post-Newtonian and 
perturbative approaches. In 1979 Smarr j2| presented his pioneer work on the numerical attack to the full Einstein 
equations. Subsequently the field developed rather slowly, in part due to a shift of the research interests toward 
other fields and in part by the lack of concrete experimental support. The situation changed dramatically in the 
early nineties with the launch of the interferometric detectors more or less simultaneously in Europe and the USA. 
Renewed efforts have been devoted to the full numerical approach in the form of the "Grand Challenge alliance" 
||, to the post-Newtonian approximation Q], and the perturbative approach. It was the perturbative method which 
in 1994 lead to the remarkable "close limit approximation" ||. Its, perhaps unexpected, success in describing the 
final stage of the binary black hole coalescence in terms of a single perturbed black hole was somewhat perturbing 
Q . A complementary approach to this approximation is to consider not two black holes of similar masses (here the 
close limit gives its optimal approximation), but instead one much less massive than the other. This, again, allows 
us to treat the problem as a perturbation around a single black hole, although here the perturbation parameter is 
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no longer the separation of the holes, but the mass ratio. This problem was recently revisited and solved for a finite 
initial separation of the holes However, both the close and the particle limits have been studied as perturbations 

around nonrotating (Schwarzschild) black holes. We know that astrophysical black holes are most likely to be rotating, 
possibly with considerable angular momentum. This means that we should take as the background metric the Kerr 
rather than the Schwarzschild geometry. 

Using the Newman-Penrose formalism, Teukolsky fic[| succeeded in giving a decoupled equation for perturbations 
around a Kerr background. Notably, this equation can be separated in all its variables. The Teukolsky equation 
has been studied and solved in the frequency domain for the gravitational radiation generated by a particle infalling 
from infinity or in circular orbit around a Kerr hole |Tl]]. Nakamura and Sasaki have given a modified version of 
the radial Teukolsky equation in the frequency domain that is better behaved (see Ref. |13| though) and suitable for 
numerical integration. Different trajectories of a particle infalling towards a Kerr black hole have been studied in this 
way [ ft4| . The disadvantage of this approach, however, is that, in general, it is much more efficient to perform the 
numerical integrations in the time domain. This was concretely our experience in Ref. g], where in order to reproduce 
in the time domain the results of Ref. jl]] obtained in the frequency domain, the running time of the program was 
reduced by three orders of magnitude. 

A successful code for the numerical integration of the Teukolsky equation in the time domain was built up only 
very recently fl. ft has been used to study quasinormal mode excitation, superradiance, and the power law behavior 
of the (late time) gravitational radiation tails on the Kerr background. Also, consistency checks similar to some 
of the material presented here, were performed during the development of the code. Given the success of the close 
limit approximation for the head-on collision of black holes, it is very interesting to study the same problem in the 
astrophysically more realistic situation of two merging black holes in an inspiral orbit. Presumably, this situation can 
be studied in terms of a single, highly rotating, perturbed Kerr black hole. To complete the solution of this problem, 
we have to provide the evolution code with consistent initial data. There are two kind of problems that we have 
to address now. Firstly, what initial data do we take? Here the problem is that, so far, all initial data available in 
the literature representing binary black holes assumed the initial three-geometry to be conformally flat. But since 
a f=constant slice of the Kerr geometry (at least in Boyer-Lindquist coordinates) is not conformally flat, we cannot 
consistently use those initial data. Only very recently a first step towards a different approach has been given by 
assuming the perturbed metric to have the Kerr-Schild form |Tq| . A second practical problem appears even if one 
had consistent initial data. How to impose them in terms of the initial Teukolsky waveform, 'J, in order to start 
the evolution? We have begun to discuss this problem in Ref. |p^| . The problem here is that VP is built up as a 
contraction of the perturbed Weyl tensor with the background (Kerr) null tetrad. This contraction, when written 
in terms of metric perturbations, depends explicitly not only on the hypersurface quantities (three- metric ^gij and 
extrinsic curvature Kij), but also on the lapse function N and the shift vector N l . We then have to rewrite '5 in 
terms of ^gij and Kij for t = (i.e. on the initial hypersurface) by making use of not only the constraints but also 
the evolution equations. In practice this task proved to be nontrivial. This is in contrast with the description, mostly 
after the work of Zerilli jl7| and of Moncrief |lS(| , of metric perturbations around a Schwarzschild hole. In this case, 
the Moncrief waveform 4>m depends explicitly on the perturbed three-geometry only and dt4>u can be written in terms 
of the extrinsic curvature only. (There are actually two waveforms representing even and odd parity perturbations 
and they correspond, respectively, to the real and imaginary parts of the Teukolsky waveform ty.) Our aim is to bring 
the Newman-Penrose- Teukolsky formulation to an equally nice footing. One first step that can be done is, at least 
in the limit of a nonrotating background, to find a relation between "J/ and 4>m on the initial hypersurface. Such a 
relation was found in Ref. jl6| making use of the simplifications introduced by choosing the Regge- Wheeler gauge and 
the possibility to restore the gauge invariance at the end of the calculation. 

When one speaks of transformations between solutions of Teukolsky and Zerilli (or Regge- Wheeler) equations, one 
immediately associates with them the name of Chandrasekhar for his contributions to the understanding of these 
relations. In the next section, we extend Chandrasekhar's transformations to the time domain in order to relate to 
</>m, and then be able to rewrite the Teukolsky waveform exclusively in terms of hypersurface quantities. We test this 
relation, that holds for every hypersurface, by numerically evolving initial data corresponding to the close limit in 
two independent ways: by the Teukolsky equation with the code detailed in Ref. [|J and by the Zerilli equation. We 
then compare waveforms, making use of the above relations, at later times, and check their agreement numerically. In 
turn, this provides further support about the correctness of the numerical code for evolving the Teukolsky equation. 
The numerical results are displayed in Sec. III. We conclude this paper with a brief discussion of our results and their 
generalization to the rotating background case. 

II. CHANDRASEKHAR TRANSFORMATION IN THE TIME DOMAIN 
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A. Teukolsky equation 



Let us briefly review the Newman-Penrose- Teukolsky description of perturbations around the Kerr metric. Gravi- 
tational perturbations with spin-weight s — ±2 are compactly written in terms of contractions of the Weyl tensor 

/9 _4 * 4 = -p^^Cnrnnfh for S = -2 

*o = ~Ci mim for s = +2 ' 



(2.1) 



where an overbar means complex conjugation, p = l/(r — ia cos 9), and we have considered Boyer-Lindquist (i, r, 9, <p) 
coordinates. This field represents either the outgoing radiative part of the perturbed Weyl tensor, (s = —2), or the 
ingoing radiative part, (s = +2). The components of the Kinnersley null tetrad Jjl| are given by 
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where M is the mass of the black hole, a its angular momentum per unit mass, S 
r 2 — 2Mr + a 2 . The source term T is built up from the energy- momentum tensor JlOj] - 



B. Zerilli equation and Moncrief waveform 



There is a historically independent way of computing gravitational perturbations around a Schwarzschild black 
hole (it can be likewise extended to any spherically symmetric background), developed mainly by Regge and Wheeler 
pof , Zerilli [0 and Moncrief jl8| . In that formalism, the two degrees of freedom of the graviton are represented by 
two scalar quantities 4>m, the even and odd parity waveforms. After decomposition of the angular part in terms of 
spherical harmonics Y[ n {9, (f>), 4>m satisfies a wave equation 



d 2 



d 2 



Ve{r)<f> u = S t (r,t) . 



(2.4) 



dt 2 dr* 2 

Here r* = r + 2M\n{r/2M — 1), Si is the contribution of the source terms, and Vt is the potential due to the curved 
background (slightly different for the even and odd parity waves). 

The even and odd parity waveforms in terms of metric perturbations in the Regge- Wheeler notation p3] take the 
form 
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(2.5) 



(2.6) 

and t is the multipole index. The field 0m explicitly depends only on the three-geometry. Likewise, one can write 
dt4>M exclusively in terms of the extrinsic curvature, using the same functional form of (j>M as above pl|: 
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C. Relation between waveforms 



The Teukolsky function can be decomposed into angular modes [ fl6| 

tt 4 (i,r,M)=£;¥i(t,r) -2Y e m (0A), (2.8) 

tin 

where we have made use of the fact that the spin- weighted spheroidal harmonics, in the case auj = 0, reduce to the 
spin-weighted spherical harmonics, i.e. e ml ^-2S" l (9,aoj = 0) = -2Y™ (8 , <fi) . 

Chandrasekhar transformations provide a differential operator (first order in frequency domain, but second order in 
time domain) that links solutions of the Teukolsky equation with solutions of the Zerilli or Regge- Wheeler equations 
and vice versa. Chandrasekhar found his transformations in the frequency domain, but in the nonrotating background 
case they can be easily extended to the time domain as (see Eqs. (3.353) and (3.345) of Ref. p2[) 

{ -Imltft, r, 0,f) } = E 7 [~ 2 ^* + 2dr.r. - V ± (r) + W ± (r) (d r « T d t )} M -^(9, <j>) , (2.9) 

where C = —\J (I — \)l{t + 1)(£ + 2)/16. This equation has to be understood as the real part of ^4 being equal to 
the even parity counterpart (labeled as "+"), and the imaginary part of — dt^4 equal to the odd parity terms (labeled 
as "— "). In the above equation V ± (r) denote the Zerilli pj and Regge- Wheeler M potentials, respectively, and 

Chandrasekhar also gives the inverse transformation for </>m in terms of a differential operator acting on ^4 (see Eqs. 
(3.319) of Ref. |§). 

We have checked analytically, in the close limit case (see below), that Eq. (2.9) leads to exactly the same relation 
as Eq. (3.6) of Ref. El for all times. 



III. NUMERICAL TESTS 

A. Close limit Approximation 

One of the more outstanding results of perturbation theory in the last years has been its application to the so called 
"close limit approximation" . The basic idea is to regard the collision of two black holes in its merger stage as a single 
perturbed black hole. When applied to two equal mass holes, this approximation gave excellent agreement with the 
full, nonperturbative, numerical computations, even going to not so small separations (5). The simplest application of 
this method considers two black holes in head-on collision. The resulting final black hole will possess no spin and can 
then be studied as a perturbation around a Schwarzschild metric. For the two black holes starting from rest at given 
separation, we have analytic expressions for the conformally flat initial geometry. In the Regge- Wheeler notation for 
the metric perturbations the only nonvanishing components of the three-metric, H% = K, are given by 

B a {t = Q t r,ff) = -^§s ± j£Z {Zo /l£)'Y?{e) . (3.1) 



2RJ l=2A 



where Zq is half the distance between the two (equal masses) holes in the conformal space. The above equation holds 
for the Brill-Lindquist data p3fl . The corresponding analysis for the Misner is equivalent to the one above replacing 
(Zq/M) £ — * Aki. The results of the Brill-Lindquist case |2j| can be compared with those of Ref. [Q, while for the 
case of Misner data |25j one can see Rcfs. H^6|. For the case of two different masses, see Ref. 



To compute </>m from Eq. (2.5) we further make the identification of the conformal radial coordinate with the 



Schwarzschild isotropic coordinate, i.e. 

R= i Vr-2M) 2 . (3.2) 

This allows us to compute th e fo rm of (here only even parity waves are generated) and from 4>m we compute 
^4 making use of the relations ( |2.9|) . Also, in this case 

1M 

d t ^ i (t = 0) = - — * 4 (t = 0). (3.3) 
We now have the explicit form of the initial data for evolving the Teukolsky equation. 
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B. Numerical results 



While gravitational perturbations can be described in terms of either sign of the spin-weight parameter s with 
\s\ = 2, it was found more convenient from the numerical point of view to deal with the s = —2 waveform (Refs. 
jl],^8|). In particular, we make a further rescaling of '5 4 and mode decomposition in the azimuthal coordinate by 
defining 

SftrV^-j-g-tfaftr*,*,*), (3.4) 



where 

1 1, . ■= il, , 4- 
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The outgoing part of <E> satisfies the asymptotic conditions 
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lim |$| ~ 1 , (3.7) 
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while ingoing solutions, i.e. those propagating towards the black hole, are characterized by 

lim |$| - 1/r 4 , (3.9) 

r* — »+oo 

lim |$| = . (3.10) 

r* — > — 00 

In the following figures we display the results obtained from runs of the Teukolsky code described in Ref. JlJ with 
grid spacings of 2 5t = Sr* — 0.05, 0.1, and 59 = ir/64, tt/32. (We have set the mass of the Sc hwa rzschild black hole 



to unity). For the comparisons with the Teukolsky evolution, the (1 + 1) Zerilli equation (Eq. (2.4) with S( = 0) was 
evolved with the same temporal and radial grid spacings as the Teukolsky equation. Previously, convergence tests 
have been performed by looking at the numerical errors for different steps of integration compared to exact solutions 
obtained by plugging a convenient analytic function into the Teukolsky equation and taking into account the resulting 
(artificial) source. In this way it is shown that the convergence is (almost) quadratic. The code can be used for any 
value of the multipole £, but in practice we used it only for £ = 2 and 4. This is so, because for the head-on collision 
of equal masses holes only even values of £ appear, and the mode £ — 2 already contributes with over 90% of the 
total radiated energy. If one considers black holes of unequal masses, one would obtain contributions from odd £, 
but still the radiation is strongly dominated by the £ = 2 mode. It should be noted that this consideration of pure 
multipoles in the sense of the spherical harmonics is only possible for nonrotating black holes: In the case of a rotating 
background, it is first of all impossible to generate initial data that represent a pure multipole specified by I and m. 
But even if we started the evolution with a pure multipole, then other modes would be generated in the course of the 
evolution. Therefore different multipoles will be present in the evolution. 

Figure 1 depicts the waveforms for £ = 2 at the initial time t — 0. We show together the initial Moncrief 
waveform 0.001 x 0m(^ = 0, r*) for the evolution of the Zerilli equation and the initial data for the Teukolsky equation 
$>(t = 0,r*,7r/2) and dt&(t = 0, r*,ir/2). (The observer is located in the equatorial plane, i.e. 9 = n/2.) Here it is 
worth to rem ark the fact that dt&(t — 0, r*,7r/2) 7^ even for time symmetric data (where dt(f>M(t = 0, r*) = 0), i.e. 



see Eq. (3.3) 



In Figure 2 we depict the initial configuration for the £ = 4 mode. 

Figure 3 shows the early behavior, as a function of the radial coordinate r* , of $(f = 3, r*,n/2) as computed in 
two different ways with Sr* — 0.05, and 59 = 7r/64: i) From the evolution of 4>wi{t — 0, r*) with the Zerilli equation 
and then using relation ( |2.9| ) to build up * 4 at the given time, and ii) by directly evolving ^(i = 0, r* , ir/2) with the 
Teukolsky equation using the code of Ref. [|l[ with the parameter a set to zero. There is an agreement between the 
two evolution methods for £ = 2 waveforms within 2% of the maximum amplitude (expected numerical error). 

Figures 4 and 5 show the same comparison of waveforms at a later time (t = 10 in units of M = 1), for £ = 2 
($r* — 0.1, 59 = 7r/32) and £ — 4 (5r* = 0.05, 59 — 7r/64), respectively, in order to check the agreement of the two 
methods of computation for different multipoles. We observe again that the curves accord within 2% of the amplitude. 
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Figure 6 shows that the agreement of the waveforms evolved with the Zerilli and Teukolsky equations continues to 
hold until late times. (Here we display curves for t = 50 obtained from an evolution with Sr* = 0.1 and 69 = ir/32). 
Numerical errors are limited to 2% of the maximum amplitude of the waves. 

These small errors have been obtained with short integration times and not so small steps in the coordinates. If we 
take them as representative of the errors in an astrophysically realistic computation we confirm both the reliability 
and efficiency of the code we used for integrating the Teukolsky equation. 

Finally, in order to have a complete description (in the context of the Teukolsky approach) of the close limit 
approach for head on collisions, we show in Figure 7 the time dependence of $ for £ = 2 as it would be seen by an 
observer located at r* = 50, 9 = ir/2. 

IV. DISCUSSION 

In this paper we reviewed the question of the imposition of Cauchy data for the Teukolsky equation. In order to 
write the Weyl scalar ^4 entering into this wave equation in terms of pure hypersurface quantities, we related ^/ to 
the Moncrief waveform 0m , that already possesses the property of depending on the three-metric and the extrinsic 
curvature only. This relation was established via a Chandrasekhar transformation extended to the time domain. 
As the Chandrasekhar transformation never makes use of any gauge fixing, the transformation between waveforms 
(already gauge invariant) and the subsequent results are gauge invariant. We checked analytically that these results 
are equivalent to those obtained in Ref . . Using this relation to impose initial data we have been able to reproduce 
the results of the close limit approximation ||, for the first time, by integrating the Teukolsky equation instead of 
the Zerilli one. In the astrophysically interesting case of the close limit, we have also been able to check numerically, 
that Chandrasekhar transformations give the same relations as found in Ref. fl6|| . Furthermore, since these relations 
hold for any hypersurface, we could test their correctness at any time by evolving the Zerilli and Teukolsky equations 
independently. We thus gained reliability about the final expressions obtained in both, this paper and in Ref. [ jl6| . 
In turn, we also gained further confidence on the numerical code developed in Ref. [jlj, that focused on the late 
time behavior of the radiation. Here the code is tested in the regime where the initial data influence is still of some 
importance as opposed to the previous studies of the gravitational tails. 

All the tests that we performed have been made for perturbations around a Schwarzschild, i.e. nonrotating, black 
hole. This allowed us the comparison with the Zerilli-Moncrief formalism. In the general case of nonnegligible rotation 
of the background hole, we do not have this counterpart to compare with and one has no other option than to evolve 
the Teukolsky equation. However, as we stated in the introduction, to the present day, we do not know explicitly 
what initial data to evolve nor how to impose them to build up "J and dt^ ■ There is perhaps an intermediate step 
we can study with the present techniques which is the case of slow rotating holes that still allow to be studied as 
perturbations of the Schwarzschild metric |^9| . 

ACKNOWLEDGMENTS 

We thank Pablo Laguna, Hans-Peter Nollert, Philippos Papadopoulos, and Jorge Pullin for helpful discussions. 
W.K. was supported by the NSF grant PHY-95-07719 and by research funds of the University of Utah. C.O.L. is a 
member of the Carrera del Investigador Cienti'fico of CONICET, Argentina and thanks FUNDACION ANTORCHAS 
for partial financial support. 



G 



W. Krivan, P. Laguna, P. Papadopoulos, and N. Andersson, Phys. Rev. D 56, 3395 (1997). 

L. L. Smarr, in Sources of Gravitational Radiation, ed. L. L. Smarr, Cambridge, Cambridge Univ. Press, pp 245-274 (1979). 

Proceedings of the November 1994 meeting of the Grand Challenge Alliance to study black hole collisions may be obtained 

by contacting E. Seidel at NCSA (unpublished). 

S. Droz and E. Poisson, Phys. Rev. D 56, 4449 (1997). 

R. H. Price and J. Pullin, Phys. Rev. Lett. 72, 3297 (1994). 

R. J. Gleiser, C. O. Nicasio, R. H. Price, and J. Pullin, Phys. Rev. Lett. 77, 4483 (1996). 
C. O. Lousto and R. H. Price, Phys. Rev. D 55, 2124 (1997). 
C. O. Lousto and R. H. Price, Phys. Rev. D 5 6, (1997). 



C. O. Lousto and R. H. Price, ^r-qc/970802g Phys. Rev. D 57, (1998). 
S. A. Teukolsky, Astrophys. J. 185, 635 (1973). 

S. L. Detweiler, Astrophys. J. 225, 687 (1978); S. L. Detweiler and E. Szedenits, Astrophys. J. 231, 211 (1979). 

M. Sasaki and T. Nakamura, Phys. Lett. 89A, 68 (1982). 

M. Campanelli and C. O. Lousto, Phys. Rev. D 56, 6363 (1997). 

T. Nakamura and M. Sasaki, Phys. Lett. 89A, 185 (1982); T. Nakamura and M. Haugan, Astrophys. J. 269, 292 (1983). 
Y. Kojima and T. Nakamura, Phys. Lett. 96A, 335 (1983); Y. Kojima and T. Nakamura, Phys. Lett. 99A, 37 (1983); Y. 
Kojima and T. Nakamura, Prog. Theor. Phys. 71, 79 (1984); Y. Kojima and T. Nakamura, Prog. Theor. Phys. 72, 495 
(1984); T. Nakamura and M. Sasaki, Gen. Rel. and Gray. 22, 1351 (1990 ). 
N. T. Bishop, R. Isaacson, M. Ma haraj, and J. W inicour, gr-qc/9711076 . 



M. Campanelli and CO. Lousto, gr-qc/9711008 



F. J. Zerilli, Phys. Rev. Lett. 24 737 (1970). We have corrected the overall sign on the Einstein tensor in Eqs. (C6a)-(C7g). 

V. Moncrief, Ann. Phys. (NY) 88, 323 (1974). 

W. Kinnersley, J. Math. Phys. 10, 1195, (1969). 

T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 

A. Abrahams and R. H. Price, Phys. Rev. D 53, 1963 (1996). 

S. Chandrasekhar, The Mathematical Theory of Black Holes, Oxford Univ. Press, New York (1983). 
D. R. Brill and R. W. Lindquist, Phys. Rev. 131, 471 (1963). 
A. Abrahams and R. H. Price, Phys. Rev. D 53, 1972 (1996). 
C. Misner, Phys. Rev. 118, 1110 (I960). 

P. Anninos, R. H. Price, J. Pullin, E. Seidel, and W.-M. Suen, Phys. Rev. D 52, 4462 (1995). 

Z. Andrade and R. H. Price, Phys. Rev. D 56, 6336 (1997). 

W. H. Press and S. A. Teukolsky, Astrophys. J. 185, 6 49 (1973). 

R. J. Gleiser, C. O. Nicasio, R. H. Price, and J. Pullin, |gr-qc/9711096. 



7 





8 



0.005 








Teukolsky evolution ♦ 
Zerilli evolution + 


0.000 






















\ 






-0.005 




& 

% 
% 


* 
* 






% 
* 
* 
* 
* 
* 
* 


* 
* 
* 




-0.010 




* 

t 
t 


* 
* 

* 








1 





-20 



-15 



-10 



10 



15 



20 




FIG. 3. Early radial dependence of <l>(i = 3,r*,n/2) for the 1 = 2 multipole computed from the evolution of cj)M.(t,r*) via 
the numerical integration of the Zerilli equation and the transformations (2.9), and from the evolution of the Teukolsky code. 
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FIG. 4. Later radial behavior of $(t = 10, r*,n/2) for the I = 2 multipole computed from the evolution of (f>m(t,r*) via the 
numerical integration of the Zerilli equation, and from evolving the Teukolsky code. 
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FIG. 5. Later radial behavior of $(t = 10, r*,n/2) for the ^ = 4 multipole computed from the evolution of (f>m(t,r*) via the 
numerical integration of the Zerilli equation, and from the Teukolsky code evolution. 
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FIG. 6. Late radial behavior of &(t = 50, r* , tt/2) for £ — 2 as computed through the integration of the Zerilli equation and 
through the Teukolsky equation. 



12 



0.04 




FIG. 7. Waveform for 1 = 2. The function $ is plotted as a function of the coordinate time t, for a detector located at 
r* = 50, 6 = tt/2. 
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